Population and Phase Coherence during the Growth of an Elongated Bose-Einstein 

Condensate 
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We study the growth of an elongated phase-fluctuating condensate from a non-equilibrium ther- 
mal cloud obtained by shock-cooling. We compare the growth of the condensate with numerical 
simulations, revealing a time delay and a reduction in the growth rate which we attribute to phase 
fluctuations. We measure the phase coherence using momentum Bragg spectroscopy, and thereby 
observe the evolution of the phase coherence as a function of time. Combining the phase coherence 
results with the numerical simulations, we suggest a simple model for the reduction of the growth 
rate based on the reduction of bosonic stimulation due to phase fluctuations and obtain improved 
agreement between theory and experiment. 
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The non-equilibrium path to Bose-Einstein condensa- 
tion is a complex process in which atoms accumulate in 
the ground state of the system and long-range phase co- 
herence develops, resulting in a strong suppression of 
density fluctuations and a uniform phase. The kinet- 
ics of condensate formation has long been a subject of 
theoretical study^giving rise to a number of conflicting 
predictions (see |l| for a review). Quantitative theories 
have been formulated to model the condensate formation 
process in a harmonic trapping potential 0, Q. How- 
ever, a limitation of these models is that the condensate 
is assumed to grow adiabatically in its phase coherent 
ground state. On the other hand, for a homogeneous 
system, Kagan et al. Q proposed the appearance of a 
quasi-condensate with strong phase fluctuations which 
die out on a time scale that increases with the size of 
the system. This homogeneous system description is also 
relevant to condensate growth in hydrodynamic clouds, 
where the trapping potential can be neglected 0. Con- 
densates in highly-elongated traps, which can often be 
treated using the local density approximation, are ex- 
ected to have properties close to the homogeneous case 
, . In addition, the axially hydrodynamic regime is eas- 
ily attainable in such traps 

Experimentally, the problem of condensate formation 
has been approached by shock-cooling 0, 0, in har- 
monic traps: starting from a thermal cloud just above the 
transition temperature, rapid removal of the most ener- 
getic atoms from the trap results in an over-saturated 
thermal cloud. Subsequent thermalization leads to the 
growth of the condensate. Measurements of the growth 
of the condensed fraction in traps significantly less elon- 
gated than ours 0, Q have obtained good quantitative 
agreement with theory 0, 0], but these experiments 
did not give access to the phase coherence of the grow- 



ing condensate. Although the two-step growth curve re- 
ported in Ref. 9], and the growth of non-equilibrium, 
phase-fluctuating condensates from hydrodynamic clouds 
m Ref. H3 support the existence of a quasi-condensate 
during the initial stage of condensate formation as pro- 
posed in Ref. Q, there exists to our knowledge no quan- 
titative experimental study of the development of phase 
coherence. 

In this Letter we present an experimental study of the 
evolution of both the condensate number and the phase 
coherence during the growth of a condensate in a highly 
elongated trap in the axially hydrodynamic regime. We 
observe that the growth is slower than that predicted by 
theory and attribute this to the presence of phase fluctu- 




FIG. 1: Condensate growth curves for initial atom numbers: 
M/10 5 = 8.6(9)., 7.2(7)A, 3.4(4)o and 2.8(5)T 20]. Each 
point corresponds to an average over three experimental re- 
alizations. The solid (dashed) lines are theoretical results of 
the model of [TH|. with (without) the correction for thermal 
phase fluctuations, as described in the text. They have addi- 
tional delay times of 10, 20, 50, and 50 ms respectively. The 
decrease in No at longer times is due to three-body losses. 
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ations which are expected to be strong in our experimen- 
tal regime, even at equilibrium Q. We then use Bragg 
spectroscopy 0, ^jj to measure the coherence length 
during the growth of the condensate. Our observations 
are compatible with the scenario of Refs. 0,|f|, where a 
non-equilibrium quasi-condensate is created at the onset 
of condensation and relaxes rapidly to equilibrium with 
shape oscillations. Finally we propose a simple modi- 
fication to the numerical simulations to phenomenologi- 
cally account for the effect of phase fluctuations on the 
condensate growth rate, and find that this improves the 
agreement between th eory and experiment. 

In our experiment |l4| . we prepare a thermal cloud 
of 87 Rb atoms in the 55*1/2 \F = {,rnp = —1) state 
in an Ioffe-Pritchard trap with final trap frequencies of 
lj± = 2ttx 655(4) Hz radially and uj z = 2ttx 6.53(1) Hz ax- 
ially. Forced radio-frequency (rf ) evaporation proceeds to 
a frequency Av r f = 140 kHz above that corresponding to 
the bottom of the combined magnetic and gravitational 
potential, giving an effective trap depth of e\ = 6 /iK, 
and the rf knife is held at this value for a time varying 
from 1 s to 12 s. This ensures thermal equilibrium, and 
allows us to control the atom number N\ in the range 2.7- 
8.6xl0 5 . The resulting thermal cloud has a temperature 
T\ of about £i/10, just above the transition temperature 
T c , which varies from 400 nK to 600 nK depending on the 
atom number. 

We next shock-cool the cloud by rapidly ramping the rf 
knife to Ai/ r f = 40 kHz in 25 ms, giving a final trap depth 
cf = 1.5 /iK. The relative truncation rate e/ef — 120 s _1 
is fast compared to the axial trap frequency, but slow 
compared to the radial trap frequency. In our elongated 
geometry this shock-cooling results in a cloud transver- 
sally at equilibrium but axially out of equilibrium. The 
cloud tends towards local thermal equilibrium with a 
temperature T < T c , in a time ~ 3r co n < 10 ms 0] 
where r co n is the collision rate at the centre of the trap 
|l6| . Since the atom cloud is in the hydrodynamic regime 
axially (uj z t co \\ « 1) 0], global equilibrium is reached 
on a time scale longer than the axial oscillation period. 

In order to study the condensate growth, the cloud is 
held in the trap for a further time t after the end of the 
shock-cooling ramp, with the trap depth held constant at 
Ef . We then switch off the trap and image the cloud after 
a 20 ms time-of- flight in order to obtain the total atom 
number TV, temperature T and condensate number Nq 
0,^3- By repeating the measurements at different times 
t for the same initial conditions, we obtain a growth curve 
for the condensate number, as shown in Fig.^for various 
initial atom numbers Ni |20j. At t ~ 20 ms (depending 
on initial conditions), the atom number has dropped by 
40% and the temperature is already below T c , yet the 
condensate does not appear until later, with a delay time 
of 20-200 ms after the fast ramp. 

We have simulated the evaporative cooling and con- 
densate growth for our experiment based on the model 



described in [llj with the additional inclusion of three- 
body loss for the thermal cloud [l^. The results are 
shown in Fig. (dashed lines). This method was in 
good quantitative agreement with experimental results 
for a less-elongated system . There are no free param- 
eters in our calculations which are based on our mea- 
sured data [2(J, however the results plotted in Fig. 1 
(dashed lines) have had a delay time of 10-50 ms added 
to better fit the experimental results. The observed 
growth rate is somewhat slower than the theoretical 
prediction. However, the simulation is not necessar- 
ily valid for such an elongated system, in the hydrody- 
namic regime, nor does it take account of phase fluctu- 
ations. Indeed, our experiment is performed in an elon- 
gated trapping geometry where temperature-dependent 
phase-fluctuations can be present even at equilibrium 
B H El El S3 For our parameters, the phase- 
coherence length = 15h 2 N /l6mkBLT at equi- 
librium is smaller than the condensate half-length L by 
a factor in the range 4-10, varying inversely with the 
condensate atom number. We will return to this point 
in more detail after discussing our experimental study of 
the phase coherence. 

We measure the coherence length of the conden- 
sate during its formation via its momentum distribu- 
tion, using 4-photon Bragg Spectroscopy as described in 
Ref. At time t after the end of the shock-cooling 

ramp, the magnetic trap is switched off and after 2 ms 
of free expansion a 2 ms Bragg pulse is applied. The 
atoms are imaged after a further 16 ms time-of-flight, 
which allows separation of the diffracted atoms. The 
diffracted fraction is measured as a function of v, the de- 
tuning between the Bragg beams which determines the 
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FIG. 2: Condensed fraction No/N and momentum width Ap 
as a function oft for an initial atom number Ni = 3.8(3) x 10 . 
(Open circles correspond to the data from Fig. Q for Ni = 
3.4(3) x 10 5 .) The decreasing momentum width indicates the 
growth of the coherence length with time. The dashed lines 
are guides to the eye. Some typical error bars are shown. 
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velocity-class diffracted, to obtain the momentum spec- 
trum of the condensate. We fit a Lorentzian function 
to the measured spectra and extract the half-width at 
half-maximum (HWHM) Av = 2ki J Ap/2Trm, where m is 
the atomic mass, &l the laser wave-vector and Ap the 
HWHM of the momentum distribution. For each spec- 
trum, further images (typically 5) are taken without the 
Bragg pulses, from which the temperature T, condensate 
number Nq, and condensate half-length L are obtained. 

The evolution of the momentum width Ap for an ini- 
tial atom number TVj = 3.8(3) x 10 5 is presented in Fig. [21 
(lower panel). The corresponding condensed fractions 
are shown in the upper panel of Fig. [21 (filled circles). 
The momentum width Ap clearly decreases as the con- 
densate fraction grows to equilibrium, indicating that, as 
expected, the coherence length grows with time. Indeed, 
for a condensate at equilibrium, we know |l3| that the 
reduction in coherence length due to thermal phase fluc- 
tuations leads to a broadening of the momentum width, 
given by (a,|2J|: 



A Pc 




0.65 



(1) 



The first term accounts for the Heisenberg-limited mo- 
mentum width due to the finite size L of the con- 
densate and the second term accounts for the pres- 
ence of thermal phase fluctuations. The numerical fac- 
tors account for integration over the 3D density pro- 
file. We also correct for the finite "instrumental width" 
of the Bragg spectrometer, by introducing as in |13| a 
Gaussian apparatus function of half-width wq = 200 
Hz. This results in a theoretical momentum width : 
A Pth = A Pcqu /2 + yj (27rm/2fc L ) 2 w 2 , + (A Pcqu /2) 2 . 
Equation ^ shows that, even if the condensate were at 
equilibrium at each instant during the growth, we would 
expect the momentum width to decrease with time, since 
both L$ and L increase with the condensate atom num- 
ber. We can therefore test whether the condensate coher- 
ence follows the density in this way by comparing each 
measured momentum width Ap with the value Ap t h cal- 
culated for a condensate at equilibrium, using the param- 
eters Nq, L and T measured for each Bragg spectrum. We 
plot the ratio Ap/Ap t h in Fig. [3 for two different initial 
atom numbers. The dashed line at Ap/Ap t h = 1 indi- 
cates the value expected for a condensate always at equi- 
librium. Systematic uncertainties of 15% on this equilib- 
rium value, mainly due to the atom number calibration 
(20%) and determination of wq (10%), are represented 
by the gray band. Unambiguously, we observe that the 
ratio Ap/Apth always lies above one (dashed line), and 
decreases in time. This indicates an excess momentum 
spread with respect to a condensate at equilibrium dur- 
ing the growth. As the condensate reaches equilibrium, 
the measurement dispersion decays and the momentum 
width tends to the predicted equilibrium value. 
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FIG. 3: Ratio of measured momentum width Ap to theo- 
retical momentum width Apth) calculated for a condensate 
at equilibrium (see text) for: (a) AT; = 3.8(3) x 10 5 and (b) 
Ni = 7.2(3) x 10 5 . (Data shown in (a) corresponds to Fig.0) 
The condensate momentum width tends to the equilibrium 
value (dashed line) at long times. Some typical statistical 
error bars are shown; the gray band indicates systematic un- 
certainties on the equilibrium value. The first vertical line 
marks the onset of condensation; the second indicates the 
time at which the condensed fraction reaches (1 — 1/e) of its 
final value. 



To interpret our results, we consider the scenario pro- 
posed by Kagan et al. [4| for condensate growth in a 
homogeneous system. The early stages of growth lead 
to a quasi-condensate with non-equilibrium, long-range 
phase fluctuations in which density fluctuations are sup- 
pressed. The phase-fluctuations then decay to produce 
the true phase-coherent condensate, with a characteris- 
tic time scale which increases with the system size 
L: Tff, oc L in the collisionless regime and oc L 2 in 
the hydrodynamic regime. Although our trapped sys- 
tem differs from the homogeneous system considered in 
Ref. 0, Svistunov Q points out that this theory can 
be applied to trapped hydrodynamic clouds, where the 
trapping potential can be neglected. In this case, the re- 
sulting quasi-condensate will be out-of-equilibrium with 
respect to a global coherent motion in the trap, thereby 
exciting a q uadrupole mode similar to that observed 
in Refs. 0|. Indeed, for a higher atom number 
N\ = 8.6 x 10 5 we directly observe quadrupole oscilla- 
tions, with an amplitude (in the trap) of 12 fim and a 
decay constant of about 250 ms. Oscillations of the axial 
condensate length broaden the momentum distribution, 
thus the excess momentum widths in Fig. 01 can be at- 
tributed to quadrupole oscillations with amplitudes of 
4/xm (JVj = 3.8 x 10 5 ) and 5.5 ^m (Ni = 7.2 x 10 5 ) [H 
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and decay constants of about 700 ms and 300 ms respec- 
tively. These values are consistent with those obtained 
for A ; = 8.6 x 10 5 , assuming an oscillation amplitude 
and decay rate which increase with the atom number. 
Therefore, apart from this decaying quadrupole mode, we 
conclude that higher-order, non-equilibrium phase fluc- 
tuations have decayed within a time shorter than 100 ms 
after the onset of condensation, in qualitative agreement 
with the predictions of Kagan et al. j2^| . 

This result allows us to interpret the data of Fig. ^ 
where the observed growth is slower than that predicted 
for a phase-coherent condensate. We suggest a simple 
picture of a ID phase-fluctuating condensate as being 
made up of a number of phase-coherent domains. As the 
growth rate includes a factor (1+ A ) due to bosonic stim- 
ulation, the growth rate of each domain will be smaller 
than that of a single large condensate, thereby reducing 
the overall growth rate. We have shown that the coher- 
ence length of our condensate approximately follows the 
equilibrium prediction, L^, which can be calculated using 
the instantaneous values of Nq, L and T measured dur- 
ing the observed growth. We therefore modify the growth 
rate in the simulation, replacing the factor 1 + No with 
1 + L^No/L |2]|. This slows the simulated growth rate, 
obtaining improved agreement with the data as shown by 
the solid lines in Fig. ^ 20] . However, the added delay 
time remains unexplained. 

In conclusion, we have observed the growth of a con- 
densate in an elongated trap and have shown, by com- 
parison with numerical simulations, that the presence 
of phase fluctuations slows the growth of the condensed 
fraction. By studying the evolution of the momentum 
width during condensate growth, we have directly ob- 
served the growth of the phase coherence with time. 
Compared with that expected for a condensate at equi- 
librium, these measurements reveal a broadening of the 
momentum distribution during growth, compatible with 
quadrupole shape oscillations. Apart from this decay- 
ing oscillation, we conclude that the condensate has 
already reached the equivalent equilibrium coherence 
length within 100 ms after the onset of condensation. Be- 
cause a complete theory for our experimental situation is 
lacking, the agreement between theory and experiment 
can only be qualitative and an extension of the model 
of Ref. Q to trapped condensates, particularly in this 
quasi-lD geometry [28[, is required. 

Using the equilibrium coherence length, we accurately 
predict the slowing of the growth rate due to phase fluc- 
tuations. However, there remains a discrepancy on the 
delay time before condensation, which may indicate non- 
equilibrium excitations at early times, as predicted by 
Ref. U A measurement of the phase coherence length 
at shorter times is needed. This is exceedingly difficult in 
our system, since at short times the condensed fraction 
is too small to obtain clear Bragg spectra. Other tech- 
niques might be used instead, such as atom laser cor- 



relation measurements |29(, combined with single-atom 
detection [30I l3ll| . 
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